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The zero temperature, or quantum, metal-superconductor phase transition is studied in disordered 
systems in dimension greater than two. A effective local field theory is developed that keeps all 
soft modes or fluctuations explicitly. A simple renormalization group analysis is used to exactly 
determine the quantum critical behavior at this transition. 



PACS numbers: 74.20.-z; 64.60Ak 

I. INTRODUCTION 

Recently there has been much interest in quantum 
phase transitions. Occurring at T — 0, these transitions 
provide new insight into the possible physical phases of 
systems at low temperature, 1 The first quantum phase 
transition studied in detail was the ferromagnetic tran- 
sition in itinerant electron system at zero temperature. 
Hertz argued in 1976 that the transition was mean- 
field- like in the physically interesting dimension d = 3. 2 
This simple mean-field description was later shown to be 
incorrectiM^ The reason for this breakdown is the exis- 
tence of soft or massless modes other than the order pa- 
rameter fluctuations. These modes were being neglected 
in Hertz's theory. In disordered systems these modes are 
diffusive, and they couple to the order parameter fluctu- 
ations and modify the critical behavior Technically, if 
these additional soft modes are integrated out, they lead 
to a long-ranged interaction and a nonlocal field theory. 
It was argued that once this effect is taken into account, 
all other fluctuation effects are suppressed by the long- 
range nature of the interactions and that the critical be- 
havior is governed by a fixed point that is Gaussian, but 
does not yield mean-field exponents. 

Similar argument was used to describe the normal 
metal-to-superconductor quantum phase transition at 
T = 0& In this case the usual finite temperature super- 
conducting phase transition is driven to zero tempera- 
ture by nonmagnetic disorder^ where the additional soft 
modes come from particle-hole excitations. Again, it was 
argued that the critical behavior found at this quantum 
phase transition 6 could be exactly determined using the 
same technique as in Refs. |^|5j 

The theory developed in Refs. 000, however, suf- 
fered from one major drawback: Since the additional 
soft modes were integrated out in order to obtain a de- 
scription entirely in terms of the order parameter fluc- 
tuations, the effective field theory that was derived was 



nonlocal- and not suitable for perturbatively calculating 
effects that depend on all of t he soft modes in the sys- 
tem. The analysis in Refs. I llol lfl was therefore 1 restricted 
to power counting arguments at tree level to show that all 
non-Gaussian terms are irrelevant in a RG sense. How- 
ever, relying entirely on tree-level power counting can 
be dangerous. In particular, logarithmic corrections can 
be easily missed. Later on logarithmic corrections were 
indeed found in the description of the quantum ferromag- 
netic transition. 9 

It is the motivation of this paper to keep all the rele- 
vant soft modes and to construct an effective local field 
theory for the metal-superconductor transition so that 
the exact behavior at this quantum phase transition can 
be determined. Unlike the quantum ferromagnetic tran- 
sition discussed above, we will show that the previous 
results for the metal-superconductor transition, though 
from nonlocal field theory treatment, are still valid. The 
inherent reason for that is explained in detail. 

This paper is organized as follo ws. In Sec. [H] we use 
methods developed in Refs. Ilfjlll to derive an effective 
local theory for disordered electron systems that explic- 
itly separates massive modes from soft ones, and keeps 
all of the latter. In Sec. IIIII we give a renormalization 
group analysis of this model. In Sec. II VI we discuss our 
results. 



II. EFFECTIVE LOCAL FIELD THEORY 

A local field theory will be given in this section to 
describe the normal metal-to-superconductor quantum 
phase transition at T = 0. All relevant soft modes will 
be contained in this field theory. We start from a gen- 
eral model of interacting electrons with quenched disor- 
der and an attractive Cooperon interaction amplitude. 
We then introduce the superconducting order parameter 
and separate massive and soft modes. After integrat- 
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ing out the massive modes, we obtain a effective iocai 
field theory that describe the coupling between the super- 
conducting fluctuations and the soft or massless diffusive 
modes. 



A. Composite field theory 

The general partition function of the interacting, dis- 
ordered electrons can be given in the form of Grassmann 
fields ip and 



D$,1>] e 



(2.1a) 



with the action S being 



f f - d 

S = - J drj cZx^ Vo-(x,r) — Vv(x,t) 

- dr H{t) . (2.1b) 
Jo 

We denote the spatial position by x, and the imagi- 
nary time by r. H(t) is the Hamiltonian in imaginary 
time representation, /3 = 1/T is the inverse temperature, 
(j = 1,2 denotes spin labels. As what we have done in 
previous papers, we integrate out the Grassmann fields 
and rewrite the theory in terms of complex-number fields 
Q and A. With the help of the following isomorphism, 
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-V'itV^t 



-i>ni>2i 



(2.2 



we exactly rewrite the partition function a: 



.10.11 



Z = I D$,i/>] e 5 ^l / D[Q]S[Q-B] 

= [ D[$,il>] S+M J D[Q]D[A] e^tAW-^l 

D[Q]D[A]e A[Q ' X] . (2.3) 



Here A is an auxiliary bosonic matrix field that plays 
the role of a Lagrange multiplier. The reason to do so 
is that the rewritten action is particularly suited for the 
separation of massive and soft modes. We then decouple 
the particle-particle spin-singlet interaction by means of 
a Hubbard- Stratonovich transformation. Denoting the 
Hubbard-Stratonovich field by the partition function 
becomes 



Z= / D[Q,A,y] e 



A[Q,A,W] 



(2.4a) 



where the action 

A[Q, A, *] = Aiis 
1 



2 Tr In [G 1 - ik 



&[Q] + 4&[Q] 

Tr (Xq) 

- / dx EE E r*s(x) r *sw 

J an r=l,2 

+iyJzr\TM\ [ d^T, E rKW 

an r—1,2 

x]Ttr [(r r ® So ) Q™ Q - m+ „(x)] • 

m 

(2.4b) 

with Tr denoting a trace over all degrees of freedom, in- 
cluding the continuous real space position, and tr a trace 
over all discrete degrees of freedom that are not summed 
over explicitly. r( c ) < is the attractive Cooperon in- 
teraction amplitude. The first three terms in Eq. I|2.4b|) 
have the following forms, 



A dis [Q] 



1 



itN f t c 



dx tr (Q(x))' 



(2.4c) 



A 



(») 



/*EH) r E E 



ni,n2,m a 



x [tr ((r r ® s ) 1 
x [tr ((r r ® s ) Q 



(x))] 

(x))] , (2.4d) 



A 



(*) _ 



Finally, 



*<£(-i) r E EE 

x [tr ((r r ® Si )Q^„ 1+m (x))] 

x [tr ((T^JiJC^W)] , (2.4c) 



G^ 1 = -d T + V 2 /2m + M , (2.4f) 



is the inverse Green operator. 

Physically, the Hubbard-Stratonovich field \& can be 
related to the superconducting, or Cooper pair, order 
parameter, \& ~ iptp. 



B. Soft modes 

Now we are ready to separate the massive and soft 
modes. Our calculations follow the same procedure in 
previous papers^ Here we just simply quote the results. 
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The effective action has the form 
A[q,*,AP,AA] = A NhaM [q} + SA[AP,AA,q} 

~ / dx EE E r^(x) r *»(x) 
•> an r=l,2 

^ an r=l,2 

x^tr (r r ® So ) [Q^_ m+n (x) 
m 

4 



(SAPS- 1 )" ^ (x) 

V / m,— m+n v J 



(2.5) 



Here if( c) = ttN^T^ /8. The matrix 5 can be treated 
as S — 1 when we neglecting the coupling between mass- 
less modes q and these massive fluctuations AP and AA, 
which is irrelevant for the purpose of the current paper. 

-4nlctM is the known action of the nonlinear sigma 
model liS 



M 



A ( nt[^FQ/4]+A ( nl[^FQ/4] 

+ ^ / ^ tr (vQ(x)" 



2# / dx tr fiQ(x) 



(2.6a) 



with A^ t from Eq. lf2~4d) . ^ from Eq. (|2~Iejl . and a 

frequency matrix with elements 

fiia = (to O so) ^12 w Bl . (2.6b) 
Here Q is subject to the following constraints, 

Q 2 (x) = 1<8>t , Q f =Q , trQ(x)=0 .(2.6c) 
and then can be write in a block matrix form as 



Q 



qq\ 



q\q 



(2.6d) 



where the matrix q has elements q nm whose frequency 
labels are restricted to n > 0, m < 0. Symmetry analysis 
with Ward identities ensures that the matrix q are mass- 
less, which are diffusive in disordered systems. The cou- 
pling constants G and H are proportional to the inverse 
conductivity, G cx 1/cr, and the specific heat coefficient, 
H oc 7 = limT^o Cy/T, respectively 14 ! 15 

5„4 contains the corrections to the nonlinear sigma 
model that were given in Ref. 0. We list explicitly the 
terms that are bilinear in the massive fluctuations AP 
and AA, but do not contain couplings between the mas- 
sive modes and q, 

5A {2] = AUAP] + J dxtr (AA(x)AP(x)) 

+ - A j rfxdytr (G(x - y) AA(y) G(y - x) AA(x)) 

(2.7) 



with A { d f s from Eq. (|2~Icjl . 

Let us first integrate out AP and AA in a Gaussian 
approximation. A additional quadratic contribution in 
terms of the order parameter field iff will obtained from 
Eq. I|2.7|) and the last term in Eq. I|2.5|l . Combining it 
with the iff 2 term in Eq. (|2.5|) yields a term 



m*] = -EEEE^w 

k ct n r— 1,2 

x [l + 2r (c) x(k,O n )] r *£(-k) ,(2.8a) 

where 

x(k,O n ) =T ^ e(mna) <5ni+n 2 ,n f nin 2 (k) , 

ni ,n 2 

(2.8b) 



2wN f t c 



(k) 



(2.8c) 



is given in terms of 

^nm(k) = ty?mn(k) 1 

with 

yw(k) = T? E Gs p(P' w «) G sp(P + k, w m ) . (2.8d) 
p 

Here G sp is the saddle-point Green function from the in- 



verse of Eq. lf2Tfj) . The Theta-function in Eq. (gjjp, 
which restricts the frequency sum to frequencies that 
both have the same sign. For small frequencies and 
wavenumbers, the calculation shows that 

A G m = - EEE E r*SCk)«ar*S(-k) , 

(2.9a) 



a n r— 1,2 



with 



U 2 = l + 0(k 2 ,fi„) 



(2.9b) 



Below we will see the wavenumber and frequency correc- 
tions indicated in Eq. I|2.9bfl are irrelevant for the critical 
behavior. 

Now we can write an effective local action includ- 
ing only soft modes and the superconducting order- 
parameter fluctuations. The action has the form of 



A eS = A G [*] + -4nl<tmM + -4c[*, q] 



(2.10a) 



Here the nonlinear sigma model part of the action, 
*4nlo-m, has been given in Eqs. I|2.6[) . and A c represents 
the coupling between iff and q, 

A c [Hf,q] = iJ*T\K(°)\ E r*"(x) 

a n r— 1,2 

x]Ttr [(r r ® S0 ) QX-m+n(x)] • 



(2.10b) 
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For the simplicity, we rewrite the coupling action as 
A c [b,q] = i^TrT\Ki c )\ J dx tr (b(x)Q(xfj . (2.10c) 



Here we define a field 

612 (x) = ^ (r r ® s ) r 6 12 (x) , 



(2.11a) 



with components 

r 6 12 (k) = S^^Sw+ntrVZOt) -(2.11b) 



Higher-order corrections can be shown be irrelevant the 
same way as in Ref. Using Eq. I|2.6d|) in Eq. (|2.1Uc|l . 
and integrating out the massive P- fluctuations, obviously 
leads to a series of terms coupling ^ and q, SS? and q 2 , 
etc. We thus obtain ^4 C [^, q] in form of a series 



Ac[y, q] = Am-q + Ap_g2 + 



(2.12a) 



The first term in this series is obtained by just replacing 
Q by q in Eq. ()2M0cJl . 

T 1/2 J dx tr (b(x) q(x)) (2.12b) 



with ci = y/ir\K^\. The next term in this expansion 
yields 

Ap_ g 2 =ic 2 Vf J dx tr (b(x) q(x) q ] (x)) (2.12c) 

with C2 = Ci/16. Terms of higher order in q in this 
expansion will turn out to be irrelevant for determining 
the critical behavior at the quantum phase transition. 



where the bare two-point q vertex has the form 



0,3 r [ 2 ) ! 3 4 ( k ) 



l,2,3p(2) (k] 
1 12,34 W 



0.3 



-<5l3^24 (k 2 + Gi?fi ni _„ 2 ) + ($1+2,3+4 

x5 aia2 5 aia3 4irTG5k c , (3.1b) 
S 13 6 2 4 (k 2 + GHn ni - n2 ) + <5i_2,3-4 
xS aia2 S aia3 4TrTGK s , (3.1c) 

S13S24 (k 2 + GHn ni _ n2 ) + 5l_ 2 ,3-4 

x6 aiaa d aiaa 47rTGK t , (3. Id) 



with K s = and K t = -irN$T®/8. Note 

that there is an additional repulsive interaction, 5k c , in 
Eq. I|3.1b() . which comes from the one-loop renormaliza- 
tion of the action. 7 Wc choose to take this effect into 
account at Gaussian order. Alternately, it would arise 
as a higher order disorder effect. For a complete discus- 
sion of this term we refer elsewhere^ Here we note that 
it is this term that drives the superconducting transi- 
tion temperature to zero, and leads to a quantum metal- 
superconductor phase transition. 

If the fermion q fields are integrated out, an effective 
action including only the superconducting order param- 
eter can be obtained. In the long wavelength and low 
frequency limit, 



k n a r— 1,2 

«2(k) + 



H 



•In 



n 



Q n \+k 2 /GH 



' "7T \n n \+^/GH 



-EEEE^(k) 

k n a r— 1,2 



t + 



\K^\ 



1 



r\h 2 In OS 

0Kc m \n n \+w/GH . 



-*S(-k) 



(3.2) 



III. RENORMALIZATION GROUP ANALYSIS 
A. Gaussian Propagators 

For the purpose of renormalization group analysis, we 
first need to determine the behaviors of Gaussian Prop- 
agators. The Gaussian or second-order action can be 
written as follows, 

A^l*,q] = -EEE E r*S(k)«a(k) r *S(-k) 



k n a r— 1,2 



?^E E E^i 2 (k);rg 34 (k); g34 (-k) 



G 



k 1,2,3,4 i.r 



BiJ*T\KV IEEE r9i2 (k) % 2 (-k), 



k 12 i.r 



(3.1a) 



Here fio is a frequency cutoff on the order of the Debye 

frequency, and t — u-2 — ^ denotes the distance from 

the mean field or Gaussian critical point. ^l < - 2 - ) [^'] is the 
Gaussian order parameter field theory that was consid- 
ered in Ref. ||J 

For the coupled field theory it is straightforward to 
compute the two-point correlation functions. For the su- 
perconducting order parameter correlations we obtain 



.*«(k),*&(p)> = a kl _p* n , m a P 



with 



M n (k) 



t + 



\KM\ 



1 

iiq 



(3.3a) 



(3.3b) 
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Similarly, we find the fermionic propagators 



mined as 



k 13 (k) ^ 34 (p)) = <5 k ,_ p Sa % *lf 2 y (k) , (3.4a) 



where 



,(2)-l 



0,3^12 34 CO ~ ^13^24^ ni-na(k) — 5l-2,3-4<^Q!ia2^aic 



and 

r ( 2 )-! fk ) 



x 2nTGK^ V ni _„ 2 (k) 2?2 _ na (k) , 

(3.4b) 

-<5l3<^24^'ni-n2 (k) + <5l+2,3+4<5a 1 a 2 <5ctlCt3 

,4^TGif( c )p ni ^ 2 (k)I? 

713—714 (k) 



l + 47iTGgWln |nn|H Sg t/Gg 



(3.4c) 



Here is the spin-singlet propagator, which in the 
limit of long wavelengths and small frequencies readsi^ 



^ i] (k) 



k 2 + G(H + K( i ))n ri 



(3.4d) 



B. Gaussian level renormalization group 

The standard momentum-shell renormalization group 
(RG) technique will be employed4£*i2'i& We use 6 as the 
RG length rescaling factor, and we rescale the wave num- 
ber and two fields straightforwardly via 

k -> k'/b , (3.5a) 
*„(Jb) -> 6( 2 -^)/ 2 *;(Jb') , (3.5b) 
q nm (k) -> 6< 2 -^/V nTn (fc') • (3.5c) 

The rescaling of imaginary time, frequency, or temper- 
ature is less straightforward. In general, there are two 
different time scales in the problem, namely, one that is 
associated with the critical order-parameter fluctuations, 
and one that is associated with the soft fermionic fluctu- 
ations. Therefore, we allow for two different dynamical 
exponents, and z q . The temperature may then get 
rescaled via 



i.—XM'pi 



T ->■ b~ Zq T' 



(3.5d) 



(3.5e) 



How these various exponents should be chosen is dis- 
cussed below. 

In the tree, or zero-loop, approximation the RG equa- 
tions for the parameters in our field theory are deter- 



(3.6a) 
(3.6b) 

(3.6c) 
(3.6d) 
(3.6e) 
(3.6f) 



Note that in giving Eqs. I|3.6e|l and (|3.6f|) . the particular 
choice of T was not yet specified because it is not obvious 
if a z q or a should be used for these terms that describe 
a coupling between q and \& fields. 

If we assume the Fermi-liquid degrees of freedom to be 
at a stable Fermi- liquid fixed point, we must choose G 
and H to be marginal, which implies 



t' 


= b 2 -^t , 


1 


b- 2 


G'H'Tq 


GHTy ' 


1 




G 7 


G ' 


H% 


= b 2 -^HT q , 


c / T ,l/2 


1 In 4_T ?* -Vq 

= Cl T x ' 2 b 2 


c' 2 T' 1 ' 2 


= C2 T 1 /V d+S ^*~ 2 '"' 



v q = 

z-q, = 2 
z„ = 2 



(3.7a) 
(3.7b) 
(3.7c) 



Here we find that two dynamical exponents, z$ and z q , 
have the same value, which is different from the ferro- 
magnetic systems^ We further choose 

m = 2 , (3.7d) 

which is implied within the logarithmic structure of Eq. 
(I3.3bl). With these choices, we find that 



/ 

c 2 



C-2 



(3.7e) 



As in the ferromagnetic case, there is a critical fixed point 
where c\ is marginal, and the fermions are diffusive, with 
exponents given by Eqs. (|3.7|) . However, in contrast to 
the magnetic case, the coupling constant C2 of the term 
is RG irrelevant for all d > 2, and so are all higher 
terms in the expansion in powers of q. We therefore con- 
clude that the Gaussian critical behavior is exacts No 
additional logarithmic corrections exist here. The most 
important technical difference that leads to the irrele- 
vance of c 2 for this quantum phase transition, while for 
the quantum ferromagnetic transition it was marginal, 
is that the time scales for the order-parameter fluctua- 
tions and the fermions, respectively, are the same^ This 
renders inoperative the mechanism that led the possibil- 
ity of C2 being marginal as in the ferromagnetic case. 
Physically, the very long range interaction between the 
order-parameter fluctuations stabilizes the Gaussian crit- 
ical behavior. This is in agreement with the fact that 
long-ranged order parameter correlations in classical sys- 
tems stabilize mean-field critical behaviori2£ 

As noted above, Eq. (|3.7e|l implies that the Gaussian 
theory gave the exact critical behavior. For complete- 
ness, the critical exponents, including logarithmic terms, 
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are critical behavior at the quantum metal-superconductor 

2 phase transition. In contrast to the disordered ferromag- 

^ — 2 — ^ n ; (3.8a) netic case studied earlier, we showed that the previous 

Ink results obtained with a nonlocal field theory were cor- 

^ _ In b red. The reason is that the two dynamical exponents, 

lnlnb 2 z<s, and z q , are exactly the same for the disordered 

7=1 . (3.8c) metal-superconductor quantum phase transition, which 

physically comes from the strong additional soft modes, 

When inoowe have 77* = 2 and v = 00. r particle-hole excitations, at zero temperature. 



IV. CONCLUSION 

We have investigated the quantum metal- 
superconductor phase transition in the present paper on 
the basis of an effective local field theory. With a simple 
rcnormalization group analysis, we have determined the 
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